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FOREWORD 


This report presents results of the expansion and improvement of the 
FORMA system for response and load analysis. The acronym FORMA stands 
for FORTRAN Matrix Analysis. The study, performed from 16 May 1975 
through 17 May 1976 was conducted by the Analytical Mechanics Department, 
Martin Marietta Corporation, Denver Division, under the contract NAS8- 
31376. The program was administered by the National Aeronautics and 
Space Administration, George C. Marshall Space Flight Center, Huntsville, 
A^labama under the direction of Dr. John R. Admire, Structural Dynamics 
Division, Systems Dynamics Laboratory. 

This report is published in seven volumes: 

Volume I “ Programming Manual, 

Volume IIA - Listings, Dense FORMA Subroutines, 

Volume IIB - Listings, Sparse FORMA Subroutines, 

Volume lie - Listings, Finite Element FORMA Subroutines, 

Volume IIIA - Explanations, Dense FORMA Subroutines, 

Volume IIIB - Explanations, Sparse FORMA Subroutines, and 
Volume me - Explanations, Finite Element FORMA Subroutines. 
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ABSTRACT 


This report presents techniques for the solution of structural 
dynamic systems on an electronic digital computer using FORMA (FOR TRAN 
^trix ^alysis) . 

FORMA is a library of subroutines coded in FORTRAN IV for the effi- 
cient solution of structural dynamics problems. These subroutines are 
in the form of building blocks that can be put together to solve a large 
variety of structural dynamics problems* The obvious advantage of the 
building block approach is that programming and checkout time are limi- 
ted to that required for putting the blocks together in the proper order 

The FORMA method has advantageous features such as: 

1* subroutines in the library have been used extensively for many 
years and as a result are well checked out and debugged; 

2. method will work on any computer with a FORTRAN IV compiler; 

3. incorporation of new subroutines is no problem; 

A. basic FORTRAN statements may be used to give extreme flexi- 
bility in writing a program. 

Two programming techniques are used in FORMA: dense and sparse. 
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SUMMARY 


The formulation and solution of most structural dynamics problems 
involves the use of matrix analysis and an electronic digital computer. 
Matrix analysis is used because it allows complicated arithmetical oper- 
ations to be formulated systematically and provides a compact form of 
bookkeeping. The electronic digital computer is used in the solution 
of the problem because of its low cost per calculation. 

After the analyst has formulated a problem in matrix notation, he 
is faced with the practical consideration of obtaining numerical an- 
swers using numerical input to the equation. The analyst must there- 
fore translate (i.e., program) the equations into a form recognizable 
by the computer. Two computer programming approaches are available to 
the analyst. One is to program the computer to solve a specific type 
problem using a basic programming language such as ALGOL or FORTRAN. 

This approach can yield a very efficient computer program but the de- 
velopment of such a program is very time consuming- Thus, such an ap- 
proach is practical only if the program will be used extensively. The 
second approach involves a library of matrix analysis operations in 
subroutine form that allows the analyst to set up his own program using 
a "building block" concept. This second approach allows the acquisition 
of quick results from problems of quite different types and is the ap- 
proach considered in this report. 

The validity of the second approach becomes evident from a study 
of structural dynamic analysis methods. This study reveals that for 
most types of problems, the mathematical operations required for solu- 
tions are limited in number. Thus, these mathematical operations can 
be programmed in the form of computer subroutines resulting in a library 
of "building blocks" that can be put together to solve a large variety 
of structural dynamics problems. The obvious advantage of the building 
block approach is that the only programming and checkout time required 
is putting the necessary blocks together in the proper order. 

The building block approach described in this report uses FORTRAN 
call statements with subroutines from a library of subroutines entitled 
FORMA (FORTRAN Matrix Analysis) . Development of subroutines in the 
FORMA library was started in 1964 by engineers in the Dynamics and Loads 
Section of Martin Marietta Corporation, Denver Division, to solve a wide 
variety of sti :ural dynamics analyses of aerospace vehicles such as 
the Titan booster and Sky lab orbiting laboratory These subroutines 
were programmed specifically for the solution of small and medium size 
structural dynamics problems of up to approximately 150 degrees of free- 
dom. Since this beginning, the FORMA library has been expanded to in- 
clude the solution of large size structural dynamics problems of up to 
approximately 6000 degrees of freedom. These subroutines for the analy- 
sis of large size structures have been used by engineers in the D>'namics 
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and Loads Section in the analysis of Viking and Space Shuttle. The 
FORMA library as included here consists of over 200 subroutines. List- 
ing and explanations of these subroutines are given in Volumes II and 

III respectively. A division is made in those two volumes for dense 

programming logic subroutines, sparse programming logic subroutines and 
finite element subroutines. 

The FORMA library includes subroutines for mass matrix calculations, 
stiffness matrix calculations, vibration modal solutions, time response 
solutions as well as the basic matrix algebra subroutines. A list of 
available subroutines is given in Appendices A, B, and C of this volume. 

The subroutines in this library have been used extensively and as 
a result are well checked out and debugged. The FORMA method has advan- 
tageous features such as: 

1. method will work on any computer with a FORTRAN IV compiler. 

It has been used on the IBM 7044, IBM 7094, GE 625/635, CDC 

6400/6500, and UNIVAC 1108 with only minor modifications; 

2* computer times are reasonable; 

3. incorporation of new subroutines is no problem; 

4. basic FORTRAN statements may be used to give extreme flexi- 
bility in writing a program; 

5. an analyst can program re“*atively complex problems with very 
little programming experience; and 

6. the method of programming is closely related to the manner 
of the mathematical formulation of the physical problem. 

In conclusion, this report expands and improves the FORMA system 
for response and loads analysis by combing existing and adding new dense, 
sparse and finite element subroutines to the FORMA library. Modifica- 
tions for MSFC requirements are included where necessary 



1 


I. INTRODUCTION 


This volmne presents the progrannning techniques and sutnmarizes the 
subroutines available in the FORMA library that will enable an analyst 
to convert his matrix equations into a computer program* It is assumed 
that the analyst has a basic knowledge of Fortran. 

Using the FORMA method, a computer program is coded using CALL 
statements for the desired subroutines from the FORMA library. Two 
programning techniques (dense and sparse) are utilized to describe the 
matrices. In dense programning, all elements of the matrix, both zero 
and non- zero, are used. The maximum size of a matrix is, thus, limited 
by the core size of the computer. For example, with an available com- 
puter core size of 50,000, the maximum square matrix size is approxi- 
mately 150 (when two matrices are used). To get around this size re- 
striction, a sparse programming technique was devised. In sparse pro- 
gramning (subroutines begin with the letter *V) only the non-zero matrix 
elements are * >ed. In the sparse technique, the matrix size is nominally 
unlimited because partitions of a matrix are stored on disk. 

A list of available subroutines is included in Appendix A (dense), 
Appendix B (sparse), and Appendix C (finite element) grouped according 
tc function (e.g., input, output, algebraic calculation, etc.). 

As with all skills, the more experienced and :^killful the analyst 
is, the ^'better'* the FORMA program he will c 'de. A “better” program is 
defined to be one that has the maximum possible matrix sizes, checks 
the input data for mistakes (\diere possible), and uses the least com- 
puter time. Probably the best means of improving FORMA skills is by 
becoming familiar with Fortran capabilities through reading of a Fortran 
coding manual. It should be emphasized, however, that any FORMA program 
will work, some programs are just “better” than others. 



II. PROGRAMMING TECHNIQUE (DENSE PROGRAMMING LOGIC) 


1- Transfer of Data 

Transfer of matrix data to and from the subroutines is made by 
subroutine arguments. Transfer of page heading data is made by a 
labeled COHNON block as explained in subroutine START. 

Input matrix data for programs using dense FORMA subroutines are 
read using Subroutine READ for real numbers (a Fortran term f?r numbers 
with a decimal point) or Subroutine READIM for integer numbers. A 
sped a l-> purpose subroutine (READO is available but is not needed for 
most programs. The only other subroutines that read input data are (a) 
Subroutine START 3 cards for (1) runs number, and user’s name, (2) title 
card 1, and (3) title card 2 ; (b) Subroutine COMENT for comment cards; 

(c) Subroutine UPDATE for tape updating data; (d) Subror.tine RBTTAB for 
data defining degrees of freedom and coordinate locations for a structural 
system. No other subroutines read input data. 

Printed output data for programs using dense FORMA subroutines are 
generally obtained by using Subroutine WRITE for real numbers or Sub- 
routine WRIUM for integer numbers. Exceptions to this are the time 
response subroutines and frequency response subroutines. Here the 
volume of calculated data is too great to be transferred out of the 
subroutine and is automatically printed in the subroutine. Other ex- 
ceptions are CKMASl, CKSTFl and RBTTAB ^ich provide specially f r- 
matted output* 

In the development of FORMA it was recognized that the matrix sizza 
and rcw dimensions could be eliminated from the subroutine arguments to 
give simpler CALL statements. However, by doing this, considerable pro- 
graming skill is then required by the analyst If in his program he wishes 
to refer to a particular element of a matrix. Considering the advantages 
and disadvantages of (1) more arguments in CALL statement but easy matrix 
element referral in main program against (.) less arguments in CALL state- 
ment but difficult matrix element referral in main program, it was decided 
to use the first approach. 
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2, Coding Procedure - Sample Problem 1 

Perhaps the best means o£ demonstrating the use of FORMA is 
through a sample problem. Assume the matrix equation 

^^*NlxN3 “ (^* ^***NlxN2 ’^*NlxN2) ^^’n2xN3 

is to be programed. Matrices |PJ, lQl> and IR) are to be input 
data to the program. The answer matrix IZ) is to be printed. The 
maximum sizes expected are N1 ^ 50, N2 *45, and N3 * 60. How- 
ever, the particular sizes of N1 , N2 , and N3 will be determined 
at run time and could be any value between 1 and the maximum size 
expected . 

The following steps are used to code the program. The program 
will be written on a sheet of coding paper to facilitate keypunch- 
ing the information to cards. A typical coding sheet with tlie 
steps listed below is shown in Figure 1. 

The names for data in a program are alphanumeric, but the 
first character must be alphabetic. A first letter of I, J; K, 

L, M, or N indicates an integer, while the rest of the alphabet 
in the first letter indicates a real number. 

Step (1) - Call Subroutine START to read three input data 
cards for (1) run number and user’s name, (2) title caro 1, and 
(3) title card 2. 

Step (2) - Write the CALL statements based on the above equa- 
tion (1) using the subroutines listed in Appendix A. This is 

shown in Figure 1 where K1 is a symbol used to designate the max- 
imum size expected for N1 . Similarly for K2,N2 and K3,N3. 

Step (J) - Write the DIMENSION statements for the matrices. 

This indicates the maximum size expected for each matrix. Note 

that an intermf;diatc matrix [PPQl * 3. IP] + iQ] is formed in Sub- 
routine AABB and must be dimensioned. The numerical values for 
Kl, K2 , and K3 are also defined. 

Step (4) - Shift back to Subroutine START by using the Fortran 
statement CO TO 1 . This procedure allows for “stacked** problems. 
The run is terminated by a STOP data card (see Subroutine START 
vrlteup) after the data of the last problem. 

iil - The end of the Fortran source deck is indl uted 
with the Fortran statement END, 




figure 1. FORHA COMPUTER PRD(»AH FOR SAMPLE PROBLEM 1 
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The input data to the sjmpio pr»Mem is also writtc'n on j tod- 
inj; lorm and is shown in Fij»ure 2. The input matrices ate assutud 
to be: 


"'2x3 - 

1 . 

2. 

3. 




4. 

5. 

t. , 



">'2x3 • 

7. 

8. 

0. 




0. 

0. 

0. . 



""3x6 - 

10. 

u. 

12. 13. 

14. 

15.' 


0. 

0. 

0. 0. 

0. 

26. 


31. 

0. 

33. 0. 

35. 

0. 


The first three cards of input data contain the following in- 
formation: 

Card 1: Run number in columns 1-6. User’s name in col- 

umns 11-28. 

Card 2 Title 1 in columns 1-72. 

Card 3: Title 2 in columns 1-72. 

The input form for each matrix is: 

First Card: Matrix name in columns 1-6, Matrix row size 

in columns 7-10 (right justified) . Matrix 
column size in columns 11-15 (right Justified). 

Middle Cards: Matrix row number in columns 1-5 (right 

Justified) of data. Matrix column number 
in columns 6-10 (right justified) of data 
in -lext field. Matrix data in four fields 
in columns 11-27, 28-44, 45-61, and 62-78. 

Last Card: Ten zeros in columns 1-10. 

The last input card is STOP in columns 1-4. 






FIGURE 2. INPUT DATA FOR SAMPLE PROBLEM 1 











3. Coding a Bor tor Program 


If the analyst is satisfied with the prottram ho has coded, 
this section can bo skipped. However, if lar>»c size nui trices 
to he used or if it is desired to check the sizes oi the input 
matrices then this section sliould be consulted. 

Equivalence - If the analyst wished to increase tlu* maximun' 
expected sizes in the program o| ire 1 to KJ = K2 * K.i - loo, 
then 50,000 core locations would he required for the matrices 
a* one. If this size requirement exceeds the capacity of the com- 
puter being used, then core locations will have to be shared be- 
tween matrices where possible. This is accomplished by using For- 
tran EQUIVALENCE. Equivalenci ng is a very sensitive operation be- 
cause it is easy to wipe out numbers of a matrix before being fin- 
ished with the matrix. Mistakes of this type will not stop the 
running of the problem and can only be noticed (hopefully!) by 
"incorrect-looking" answers. 

There are several methods of equiva fencing. In the first meth- 
od, the various matrices are equlvalenced to locations in a large 
dummy matrix. In the second method, the various matrices are equiv- 
alenced to each other. The third method is a "manual equivalenc- 
ing** procedure and is recommended over the other two methods be- 
cause it is easier to code and understand. In this manual equiv- 
alenclng method, only two or three matrices are dimensioned [e.g,, 
DXMENSICN A(100,100), B(10Q,100), C(100,100)J. The entire program 
Is coded using only the nemes A, B, or C. The particular meaning 
of A, B, or C should then be given In Columns 73 thru 80. Bv 
this third method, the EQUIVALENCE statements are kept to a mini- 
mum and may not be needed at all. 

It Is advisable to equivalence only the larger matrices of a 
program. The possibility of mistake introduced by equivalenclng 
scalars and the S'-.ialler matrices is not worth the small amount of 
core that will be saved. 

To demonstrate the manual equivalencing method, assume the di- 
mensions of the program of Figure 1 are to be increased to K1 
K2 * K3 « 100. Assume that the resulting 50,0C0 core locations 
exceeds the capacity of the computer being used. 

In the first example of manual equivalencing, Subroutine MULT 
is retained In the program. Thus a minimum of three matrices will 
be needed. The resulting program is shown in Figure i. The core 
requlreioents for this program art- only 30,000 for the matrices. 

The input data of Figure 2 is still the same. The samt* results 
as the program shown in Figure 1 will still Ik- obtained. 





FORMA 













In the soronci cxamplt.* of i-.kinu-tl vqui v.i 1 \*nvi *>>•. , thu prnjir.i'’: oS 
Fi>;urc J is inocMfic-O lo u^.»■ ^mi!- * -s' :ujtriics. iii s riquiri> 

Sul»roiitiiH» MULT bo repfacoO with fiihv-r Subrout i no MULTA t>r Ml'I.'iB. 
Tho resulting program is shown in Figure 4. Tho tore roquln’**H-nf s 
lor this program aro only 20,000 for the matrKos, Tho input Ju^u 
o! Figure 2 is still the same. The same results as the program 
sliown in Figure 1 will still be obtained. 

Size Check - Any of the three programs just coded '..ill rtin 
even if there is a mistake in the matrix sizes in the input uata. 
The modi fications to the program of Figurr* 4 to cheek tho si:-.es 
of the input matrix data are shown in Figure 5. This is tlu *'hest" 
program for the sample problem given by equaticn (1), The input 
data of Figure 2 is still the same. The sa:ne results as the pro-- 
gram shown in Figure I will still be obtained. 


ORIGINAL PAGE IS 
POOR QUALErX 



FORMA 




FOMA 


1 



FIGURE 5. FORMA COMPUTER PROGRAM FOR SAMPUB PROBLEM 1 
TWO MATRICES DIMENSIONAL 
SIZE CHECKS INCORPORATED 
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Sample Problem I 

To £urtl*er illustrate the use of FORMA, a second sample prob- 
lem is coded in this section. The coding techniques described in 
Section 4 to obtain a ‘‘better** program will be used. 

In this problem, the “free-irec** mode shapes and frequencies 
of the beam shown in Figure 6 arc to be calculated and printed. 

Two degrees ol freedom, translation and rotation, are assumed at 
each of the five panel points (also called collocation points) . 

The input data to the computer program are the beam panel points, 
the beam weight distribution, and the beam stiffness distribution. 
For this sample problem the distributed rotary inertia, any con- 
cencrated weights, and the shear stiffness are ignored. 

The following steps are used to code the computer program. 

The program will be written on a sheet of coding paper to facili- 
tate keypunching the information to cards. A typical coding sheet 
with the steps listed below Is shown in Figure 7. 

As mentioned previously, the names for data in a program are 
alphanumeric, but the first character must be alphabetic. A first 
letter of I, J, K, L, M, or N indicates an integer, while the rest 
of the alphabet in the first letter indicates a real number. 

Step (1) - Call Subroutine START to read three input data 
cards for (1) run number and user's name, (2) title card 1, and 
(3) title card 2. 

Step (2) - Write the CALL statements to read in the panel 
points, weight distribution, and stiffness distribution. Checks 
on the column size are made. Write the CALL statements to calcu- 
late and write the mass and stiffness matrices, and to calculate 
and write the mode shapes and frequencies. K1 is a symbol used to 
designate the maximum number of degrees of freedom allowed. K2 is 
a symbol used to designate the maximum number of panel points al- 
lowed. K3 is a symbol used t(~ designate the maximum number of 
rows of distributed data. 

Step (3) - Write the DIMENSION statements for the matrices. 
This indicates the maximum size expected for each matrix. Even 
though this sample problem has five panel points, the computer 
program is written assuming that there could be as many as 30 pan- 
el points and thus 100 degrees of freedom. Also, a maximum of 
40 rows of distributed data is allowed by the dimension given to 
the matrix D. The corresponding values for K1 , K2 , and K) are 
defined . 
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1 2 3 4 5 

■ 

0 30 50 100 120 

Station (In.) 

(a) Panel Point Arrangtnent 



Station (In.) 

(b) Weight Distribution 



u 


Station (In.) 

<c) Stiffneaa Distribution 


Fifura 6. Bsaai for Sanpla froblen 2 
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Step (-<») - Shift hjck io Sn)>rc'Mtine START by using the Kiutran 
statemenl 00 TO 1. This procedure allows for ’’stacked** problems. 
The run is terminated by a STOP data card (see Subroutine START 
writeup) after the data of the last problem. 

Step (5) - The end of tlie Fortran source di‘ck is indicated 
with the Fortran statement END. 


The input data to sample problem 2 is also written on a cod- 
ing form as shown in Figure 8. The first three cards of input 
data contain the following information: 

Card 1: Run number in columns 1-6. User’s name in col- 

umns 11-28. 

Card 2: Title 1 in columns 1-72. 

Card J: Title 2 in columns 1-72. 


The input form for each ntatrix is: 


First Card: 


Matrix name in columns 1-6, Matrix row size 
in columns 7-10 (right justified). Matrix 
column size in columiis 11-15 (right justified). 


Middle Cards: Matri c row number in columns 1-5 (right 

justified) of data. Matrix column number 
in columns 6-10 (right justified) of data 
i.n next field. Matrix data in four fields 
in columns 11-27, 28-44, 45-61, and 61-78. 


Last Card: Ten zeros In columns 1-10. 


The matrix data consists of: 

1) Matrix of panel point stations from Figure 6(a) . 

2) Matrix of end point coordinates of the line segments 
representing the distributed weight from Fig\ire (>(b) . 

3) Matrix of end point coordinates of the line segments 
representing the distributed bending stiffness from 
Figure 6(c) . 
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The end point coor'dinates of each nonvertical straipln 1io«. lv>r 
Che distributed data is given as a row in the matrix oi di.sirih- 
uted data. Each row has the form: 

Matrix column 1 - station x^, i,e., the abscissa of tl>c 
line segment originating point. 

Matrix column 2 - station ^ , l.e., tlic abscissa of the 

line segment terminating point. 

Matrix column J - value at x^(+) , i.e., the ordinate ol the 
line segment originating point. 

Matrix column A - value at l.e,, the ordinate of the 

line segment terminating point. 

The last input card is STOP in card columns 1-4. 



ORIGINAL PAGE 
OF POOR QLALn 



FICJRS 8. INPUT DATA FOR SAMPLE PROBLOl 2 
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III. PROGRAMMING TECHNIQUE (SPARSE PROGRAMMING LOGIC ) 


1 . Transfer of Data 


Matrix data for the sparsi ibrou tines is stored on disk with a disk 

number (representing a matrix) being transferred to and from the subrou- 
tines by argument. Transfer of page heading data is made by a labeled 
COMMON block as explained in subroutine START. 

Input matrix data is read using subroutine YREAD and printed output 
data s obtained using subroutine YWRITE. 

2. Coalng Procedure 

The same example will be used here that was used for sample problem 1 
for the dense programming logic. The example is repeated from page 3: 

H Nl3dS3 “ ( ^ H NlxN2 [^] NlxN2 ) [^] N2xN3 

As before, matrices [p] , [q] , and [r] are to be input data to the 
program. The answer [z] is to be printed. 

The following steps are used to code the program. The program is 
written on a sheet of coding paper to facilitate ke> punching the in- 
formation to cards. A typical coding sheet with the steps listed below 
is shown in Figure 9. 

Step (1) - Dimension workspaces V and LV at least 3 times the largest 
row or column size expected. The larger the dimension size the faster 
the computer time. Indicate the dimension size with 

Step (2) - Set the tape names to numbers. 

Step (3) - Call subroutine START to read three input data cards for 
(1) run number and user’s name, (2) title card 1, and (3) title card 2. 

Step (4) - Write the CALL statements based on the above 2 quation (1) 
using the subroutines listed a Appendix B. 

Step^lSi - Shift back to subroutine START by using the Fortran statement 
GO TO 1. This procedure allows for "stacked" problems. The run is ter- 
minated by a STOP data card (see subroutine START writeup) after the 
data of the last problem. 

Step (6) - The end of the Fortran source deck is ividicated with the 
Fortran statement END. 
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The input data for this sparse program is identical to the input 
data previously s own on pages 5 and 6 (Figure 2) for the dense program. 

The techniques of pages 7 through 11 for coding a better dense pro- 
gram are not pertinent for a sparse program because equivalence is not 
needed. The size checks could be made with the sparse program but is 
not shown here. 
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FIGUR£ 9. FORMA COMPUTER PROGRAM, SPARSE PROGRAMMING LOGIC 





A-1 


APPFNOIX a. summary TF calling INSTPUCTIPNS - 
01 NSC FORMA SUFROUTINFS. 


IN 7HF ARGUMENTS CF THF SUGROOTINES PFLPW IT IS ASSUMFU THAT 
THFRF IS CrPFESPrNDFNCE IN SIZF ANO ROW DIMF^'SION OF 
C0MPATI*'LF MATRICFS, for INSTANCE IN SUPRaJTINE MULIt 
NRAsNRZ* N(A*NRBt NCB*NCZ* KA^KZ 

A 6H AFrUMFNT MAY ALSO FE A VAR I APIS READ WITH AN A6 FORMAT CP 
OPTAINFD WITH A DATA STATEMENT. 


LIST OF SYMBOLS ♦♦♦ 

A = INPUT MATRIX 

e = INPia MATRIX 

ALPHA s INPUT SCALAR 

BETA = INPUT SCALAR 

AVFC = INPUT VECTOR tPOW OR COLUMN MATRIX I 

IVFC = INPUT VECTOR fROW CP COLUMN MATRIX), INTf&fR 

TAP = INPl»T TAPLE IMATPIX WITH INCOMPLETE COH>MHS IN SOME ROWS) 

Z = RESULT MATRIX 

ZVEC = RESULT VECTOR 

N = SIZE (FOP VECTOR OR SCUARE MATRIX) 

NR = NUMEER OF ROWS 

NC = MiMPER OF COLUMNS 

K = ROW DIMENSION 

KR = ROW DIMENSION 

KC = COLUMN DIMENSION 

V - VECTOR WORK SPACE 

LV s VECTOR KfRK SPACE, INTECEP 

KV = V,LV DIMENSION 

NUTI = LOGICAL NUMBER OF ITH UTILITY TAPF 
NTAPF = SYMBOLIC NUMBER OF TAPF. FOR FXAMPLE , 1 


A SINCLY DIMENSIONED VARIABLE IS REFERRED TO AS A VECTOR IN THIS 
REPORT. 

A PCUBLY DIMENSIONED VARIABLE IS REFERRED TO AS A MATRIX IN THIS 
REPORT. 

A VECTOR MAY PE HANDLED AS EITHER A ROW PR COLUMN MATRIX. 

THE ROW DIMENSION OE A VECTOR IN THF APGUM^^ ‘ S OF IHE CALL 
STATEMENTS IS ANALOGOt'S Tl THE ROWS OE ThE VECTOR. THAT IS, 

IF THE VECTOR, IS HANPLFC AS A ROW, THEN KR=1 

IF THE VECTOR IS HANOLFP AS A COL, THEN KR=UIMENS ION SIZE 




A-? 


.01 MISCFLLANFDUS 

.01.01 POC'GRAH INITIALIZATION 

call start 

.Cl. 02 PROGRAM PAGF HEADING 

CALL PACFhP 

.01.03 PROGRAM PCMPOUT 

CALL ZZPOMB l6HSUGNAM,NERROR) 

.01.0^ PROGRAM COMMENTS 

CALL COMENT 

.01.05 CONVERSION 

CALL YDTOS ( A ,NUTZ , NR ,NC ,KR ,KC . V,L V,K V »NUT 11 
OSES YIN ,YINI ,Y0L!T,Y0L)T1,YPART,Z2F0MP. 

.01.06 matrix element COMPARISON 

CALL COMPAR ( A ,R FF ,N» ,NC tNI'IG tCTOL . 6HAN AME-,6HRFFNAM, 
KAfKREF) 

.01.07 TIME CHECK 

CALL TIMCHK <6HNAMCHK| 

.01. OP ORPFRING 

CALL XICRD <V,LV»LAStNN2A) 

CALL ORDALP ( IMAT.NR ,NC ,NC AL* IWM AT, KRI . KCW ) 

.01.09 MERGE NAME AND NUMEEP (FUNCTION) 

NAME (6FNAMF ,NOM ) 
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.0? INPUT 

•02. Cl REAL NUMPfP? 

CALL PLAT (A,ANP,*NCtKP,KC) 

uses INT APF ,LTAPF ,PAr,CHD ,PTAPF ,WP ITF tWT APF ,Z 2FDWF. 
♦NP,NC WILL PE PEFINFP PV INPUT. USE SYMPCLS. 

CARD INPUT 

FIPST CARD “ 2NAMF»NP,NC WITH A6tI4,If» FORMAT. 

RFMAPKS IN rCLl'MNS 16-60. 

S IN CCL 7? FCP WP1?E TAPF /NlTlALIZF, 
PLANK, REWIND, L 1ST, OR TAPFID <FC'R WRITE 
TAPE, IN CriUMNS 73-7F. 

NWTAPE IN COLUMNS 7O-R0. 

MIDDLE CARDS - DATA WITH 2IF,4F17 FOPHAT. 

LAST card - TEN ZERfFS IN CCLUI-.is- 1-10. 

TAPF INPUT 

ONF CARD - ZNAME,0 OP -LCCAT ION ,NRTAPE (W ITH - FOR 

NO WPITf 0UT),2RUNN0 WITH 
A6,I4,1«;,A6 FORMAT. 

PLANK, REWIND, OR LIST (FOP READ TAPE) 

IN COLUMNS ??-77. 

REMARKS IN COLUMNS 2E-6P. 

$ IN COL 7? FOR WRITE TA^f INITIALIZE. 
FLANK, REWINOtLIST, CP TAPEID (FOR WRITE 
TAPE) IN COLUMNS 73-78. 

NWTAPE IN COLUMNS 70-FO. 

.C2.C? INTECFR NUMPFPS 

CALL PFATIM ( IZ ,+NH,*Nr ,KR ,KC ) 

USES INTAPE ,l TAPE,PAGFF'n,PTAPF,WF ITIM,WTAPE ,2 2BOHE . 
*NR,NC WILL BE DEFINED EY INPUT. USE SYMBOLS. 

CARD INPUT 

FIRST CARD - I2NAMF,NR,NC WITH A6,14,lf FORMAT. 
REMARKS IN COLUMNS 16-6S. 

* IN CPI 72 FOR WRITE TAPE INITIALIZE. 
FLANK, REWIND, LIST, OR TAPFID ( Ff‘P WRITE 
TAPE) TN COLUMNS 73-7F. 

MWTARF IN COLUMNS 7P-FC. 

MIDDLE CARDS - DATA WITH 216,1416 FORMAT. 

LAST CARD - TEN 2FR0FS IN COLl*MNS 1-lC. 

tape INPUT 

(SAME AS SUPPOUTINf READ ALOVE). 

.02.03 OCTAL NUMFFPS 

(ALl RE ADO (Z,*NR,*NC,KR,KC) 

»NP,NC WILL BE DEFINED PV INPUT. USE SYMEriS. 

FIRST CARD - ZNAME ,N^ ,NC WITH A6,14,JS Ff>RMAT. 

REMARKS IN COLUMNS 16-6«>. 

MIDDLE CARE'S - DATA WITH 2 1 6 ,? (3X ,02C ) FORMAT. 

LAST CARD - TEN ZFROFS IN COLUMNS 1-10. 





.02 INPUT ICONTINUED) 

.02.04 alpha-numeric characters 

CALL RFADAN <12 ♦♦NRt*NC ,KP ,KC ) 

USES IMAPE»LTAPE,PACFK',RTAPr,KRlTAN,WTAPF,Z2BCME. 
♦NR,KC WILL BE DEFINED PY INPUT. USE SYMBOLS. 

CARD INPUT 

FIRST CARD - I2NAME»NR»NC WITH A6, 14,15 FORMAT. 
REMARKS IN COLUMNS 16-69. 

S IN CCL 72 FOR WRITE TaPF INITIALIZE. 

blank,rfwind,list,cr tapeid (fcr write 

TAPF) IN COLUMNS 73-7B. 

NWTAPF in COLUMNS 79-80. 

MIDDLE CARDS - DATA WITH 2I5,I0A6 FORMAT. 

LAST CARD - TEN ZEROES IN COLUWS 1-10. 

TAPE INPUT 

(SAME AS SUrPOUTINF READ AFCVE). 
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.03 OUTPUT 

• 03.01 PPINT RFAL NUMPF<?S 

CALL WRITF < A,K‘R t^C ,6HANAME-,K ) 

.03.0? PRIM IMfCEP NUMPEff 

CALL WPITIM (lAfNR tVC.6HlAN*AMF,K 1 

.03.03 PRINT ALPHA-NUMERIC CHARACTFP< 

CALL WRITAN ( 1 A *NR t NC ,AHI AN AMF ,K > 

.03.04 PLOT 

CALL PLCTl <XVEC,YMAT,NPY,NCY,XfTAPT,XrELTA,6HXNAME-t 
YNAME»PTITLF,TFSAME,TFCURV,IFLIF1, K> 

USES PLOTSS. 

HAVE -CALL PPLT (C,?HLCI- IN MP USED ONCE/PUN. 

MAX NCY = 3 

CALL PLOT? IXVEC,YMAT,NRY,NCY*6HXNAME-,6hYNAME-, 
PTTTLE,IPLnT,IYS, KY) 

MAX NCY=10 

CALL PLPT3 ICLOf ,MLUC,CCf LOC^VPLCC .PANGLF.CANGLE, 
fEDIST,IFJNUM»LRFYF ,NV IFW , IFF A t PTITLF , 

NC fNM,KC,KMI 
USES VCROSS,VrCT. 

CALL PLCTSS I YMAX, YM IN, YTOP ,TFCT » 

.03.05 PUNCH 

CALL PUNCAN C 1 A,N» ,NC,6H1 ANAMF ,K ) 
call PUNCH (a,NR,NC ,6HANAMF~,K» 

CALL PUNCHO < A,NR ,NC ,6HANAMF-,K ) 

CALL PUNCIM <IA,NF ,NC,tHlANAMF,K ) 
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.04 TAPF CPFRATIPNS 

RFWIND NTAPF AT PEGINNINC RF MAIN PPPCPAM. 

.04.0! IMITIALlZATirW 

CALL INTAPE ( NTAPF ,6HTAPF IP ! 

.04.02 REAPING 

CALL PTAPF (6HZRUNNP,6HZNAMF-.2,*NR,*Nr ,KP,KC,NTAP£ ) 
♦NP*NC WILL BE OFFINFO PY TAPE DATA. USE SYMBGLS. 
WSFS LTAPE. 

.04,03 WRITING 

CALL WTAPE ( A .NR .N*" .6HAN AMF-.K , NTA PE ) 

.04.0^ LISTING DF HE; HERS 

CALL LTAPE (NTAPF) 

.04.05 UPDATING 

CALL UPDATE 

.04.06 CCPE/TAPE DATA TRANSFER 

CALI IN (NTAPEtZtN) 

CALL PUT (NTAPF .A .N) 

CALL RWND (NTAPE ) 

CALL SKPP (NTAPE .NKEC) 
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C5 

GFNFPATICN 



05.01 

BASIC 




CALL 

ONES 

(Z,NRfNC,K) 


CALL 

SIGNA 

(ZtN«K) 


CALL 

UNITY 

(Z«N«K) 


CALL 

ZFPO 

(ZyNRyNC «K) 
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.C6 algebra 

.06.01 ECALAE PPDPUCT 

CALL ALRHAA ( ALPHA ,A ,Z ,NR ,NC ,K A2 » 

C#LL PA (P,A,2,NR,NC,KA,KZ) 

.06.02 ADDITION, SUBTRACTION 

CALL AAFP (ALPHA,A,PETA,B,Z,NR ,NC ,KARZI 
CALL PACP IP,A,0,P ,Z,NR ,NC,KA,KP,KZ» 

.C6.C3A NULTIPLICATION - GENERAL 

CALL MULT ( A,B ,Z ,NRA ,NPB,NCB ,KA2 ,KB) 

CALL MULTA (A2,P,NRA,NRP ,NCP ,KAZ,KB 1 
MAX NRB = BOO 

CALL MULTB < A,P2,NP A,NPE ,NCB ,KA ,KP Z ) 

MAX NRE - fOO 

.C6.03P MULTIPLICATION - SPECIAL 

CALL API (A,P ,2,NPA,NrA,NCP,KA,KE,K2 I 
MAX NCA = 500 

CALL AP? (A,B,Z,NPA,NCA,NCFs,KA,KE,KZ » 

MAX NCA = 500 

CALL APCl (A,B,C,2,NRA,NCA,NCB,KA,KB,KC,K2) 

MAX NCA = 500 

CALL ABC2 < A,E ,C ,2 ,NPA ,NCA,NCB ,K A ,KP ,KC ,K2 I 
MAX NCA = 500 

CALL ATPl (A,P,Z,NRA,NCA,NCP,KA,KP,KZ » 

MAX NPA = 500 

CALL ATP2 (A,P,2,NPA,NCA,NCe,KA,KB,KZ) 

MAX NPA = 500 

CALL ATPCl (A,e ,C,Z,NPA,MCA,NCP ,KA,KP,RC,KZ ) 

MAX NPA = 500 

CALL ATPC? fA,P,C,2,NRA,NCA,NCP,KA,KP,KC,K2) 

MAX NRA = 500 

CALL ATXPAl CA2 ,P ,NRP,NCP ,KA2 ,KP ) 

MAX fCP = 50C 

CALL ATXPR ( A ,P2,NR AT,NRF ,NC F ,K A,K P2 ) 

MAX NPAT = 500 

CALL ATXFPl ( A,E 2 ,NP E,NCF ,5 A,KP 2 ) 

MAX NRB = 5CC 

CALL ATXPP? (A,F2,NRP,NCE,KA,KB2» 

MAX NCR = 500 

CALL AXPAl f A2,P,NRA,NC A,KA2,KB I 
MAX NCA = 500 

CALL AXPA2 ( A2 , B ,N , K A2 ,K E ) 

MAX N = 500 

CALL AXPA3 C A2 ,P ,NRP ,NCE ,K A7 ,K E ) 

MAX NPP = ^CC 

CALL OBI (D,B,2,NRP,NCP,KP,K2 I 
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.06 ALrEFPA (CONTINUED) 

.06.04A TRIPLE MATRIX PPOnUfT - GFNFFAL 

CALL PAPT lA,F,2fNRP,NCP,KA,KPZ) 

MAX NCP = fCO (SIZE OF A) 

CALL BAFT A ( AZ»F tNPP.NCF ,KAZ ,KF I 
MAX NCP = 5C0 ISIZF TF A) 

CALL RTAP (AtP,Z,NRP,NCP,KAP,KZ) 
max NPF = 500 (SIZE OF A I 
CALL PTAPA lAZtf ,NRB,NCP.K.A2,KP) 

MAX NPB,NClr = 50C (SIZF CF MATRICES A,2I 

.C6.04P TRIPLE MATRIX PPrOUCT - SPFClAt 

CALL ETAPl (A,P,Z»NPP,NCF,KA,KP,KZ) 

MAX NRP = 500 (SIZF OF MATRIX A) 

CALL BTAPA? (AZ,P,N,KA» 

MAX M = 5CC (SIZE CF MATRICES A,P,ZI 
CALL BTAPCI (A,P,C,Z,NRP,NCE,KA,KP,KC,KZI 
MAX NRP = 50C 

CALL BTDBI (D,&,Z,MPB,NC'’,KRPfKZ ) 

MAX NRP = 5C0 

CALL PTDBCI (D,E ,C,Z ,NRP,NCF ,KR ,KC,KZ) 

CALL UTAUl (AtU,Z,N,KA,KlJ,KZ) 

MAX N 500 

CALL OTAUCl ( A ,U *C tZ ,N,K A ,K U*KC » KZ ) 

MAX N = 5C0 

.06.05 INVERSION 

CALL INVl (A,2fN,KAZ) 

MAX N = 250 

CALL INV? (AfZ,N,KAZI 
MAX N = 250 

CALL INV3 (A,ZtN,KAZ) 

USES DC0M1,INV4. 

MAX N = 250 

CALL INV4 (AtZtNtKAZI 

.06.06 SIMUITANFOUS ECUATKZNS 

CALL SMEOl (♦A**PVFC,ZVEC ,N,KA ) 

*A,PVFC ARE DFSTROYEO. 

.06.07 ETCENVALUF, EICENVFCTCR 

CALL EICA^l (♦A.ZVAl ♦ZVFC»N*FCn,KAZ I 
*A IS DESTPCYED, 

CALL EIGNIA ( ♦ A, ZVAL »Z VFC ,N1N ,( TV,K AZ ) 

♦A TS DESTPOYED. 
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,06 ALGFFRA (CONTINUED) 

.06.08 DFCCMPOSITION 

CALL nro^l (A,Z,N,KAZ) 

.06.0*» ROW, COLUMN CPFRATICt*S 

CALL ROWMLT ( AVFC,P,Z,NR ,Nf ,KFZ I 
CALL COLMLT ( AVEC , P , Z ,NR ,NC ,KPZ ) 

.06.10 VFCTOP OPFPATICINS 

CALL VDOT (VA,Ve,PROnCT,VAMAf ,VFMAG,COSAF ) 
CALL VCRCFS ( VA,VP ,VZ , VAM AG ,VPM AC,VZMAG ,S INA8 ) 
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.07 MODIFIfATICN 
.07. 01 ASSEMBLY 

TALL ASSFM ( A , 1 2 , J2 t »2 ,NFA,N'f A ,NPZ ,Nf ^ ,KA ,KZ J 

♦BE SURE Z IS DEEINED. <Ef CALL ZKl TP CLEAE Zl. 

.07.02 DISASSEMBLY 

CALL DISA U»lA,JA,2.NRA,NCA»NRZf NCZ *KA.K2 ) 

.07.03 REVISION 

CALL REVADf) C ALPhA t A tl VFC t JVfc , *2 »NP A,NCA*NRZ ,NC2 t KA » K2 ) 
♦EE SURE 2 IS DEFINED. (EG CALL ZFPP 10 CLEAR 2). 

CALL REVIJ (AZf IVfC,JVFC,NRA,NCA,NPZ,NC2,KRAZ) 

.07.0A TRANSFOSF 

CALL ATI (A,2 ,NPA,NCA,KA,KZ ) 

MAX NCA = BOO 

CALL TRANS ( A »Z tNRA f NCA , •< A ,K Z I 

.0'».CB SYMMETRIZE 

CAL'- SYMLH (AZ,N,K) 

CALL SYMUH (AZ,N,K) 

.07.06 NULLIFY 

CALL ZFPOLh (AZ,N,K) 

CALL 2ER0UH (AZfN,K) 

.07.07 DIAGONALIZE 

CALL DIAG (AVECtZtN»K2 ) 



OP INffPPOLATIPN, 01 FFFPPN'flATK N 

08.01 INTFRPCLATICN 

TALL TFRPl ( XA,X2 ,A »2 .F'XAtNXZ ,NCa, KA,K 2 ) 
CALL TFFP2 « XA ♦ XZ ,A »2 »NX A ,NX 7 ,Nf A, K A ,K 2 ) 

08.0? OIFFFRFMIATION 

CALL OIFFl (XA,XZ,A»Z,r'XA,NXZ,r^CA.KA,K2) 
CALL PIFF2 (XA,X2tA,ZtNXAtA'X2,NCA»KA,K2l 
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.CO AIRlCAr 
•09.01 LATERAL 

CALL ALCDl fPP.DI ST ,CONC,CrNVP T ,2VfC *NPP,ND,NC .KO.KC I 
.09.0? AXIAL 

CALL ALr02 (PP.OIST.CONC ♦CnNV9T,2VEC,NPP,ND,NC,KDf KC I 
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.10 

f*AS5 


.10.01 

pro 

CALL rASSl 

|PP,DMASS,0RlN,C0NC,CCNVPT,2.NPP,NDH,MritNC, 

KPM.Kni.KCtKZ) 

.10.02 

FF AM 

CALL MASf2 «PP ,DMASf tOP IN.CCNC ♦ CCNV^T, 2 ,NPP,»«OM,NOI ,NC, 
♦NZ ,KDM,KDI,KC,K2I 

♦MZ=2NPP WILL PF DFFINFD PY MAS52. USE SYMFCL. 

.10.03 

FLUlf' 

CALL MASS2A 
♦N2=2NPP+1 

|PP,OMASS,rCSM,FLFVFL,CCf!VRT,2f^PP,MCM,*t>v2t 

K0K,K21 

WILL BF DFFINFP PY MAS52A. USE SYMPCL. 
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.11 STIFF»^>ES5 


.11.01 Fcn 

CALL STIFl IPP*OAE,2tNFPtMrAF,KOAE,»c2l 


. 11 . 0 ? 


PF AP 

CALL STIF? 


♦N2=2NPP 


(PP*DKAG,nFI ,2,^'PP,^^vKAC,tNDf 1**N2*KCKAG» 
KHFI.KZ ) 

WILL FF DFFINFP EY <^T1F?. USE SYKECL. 


.11.03 


PFPUCTIPN 

CALL SPfni 
CALL SFED2 
CALL SFFC** 


1A,P,T,N.NR,1FT,KART1 
(A,P ,1 ,N,NF ,IFT,KAPT1 
(A,IV,FtTtKtNP»IFT,KAFl ) 
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.12 MPHAl 
.12. Cl JACfPl 

CALL MP['E1 (♦AMA5S,**5i1IC,W?,W,FPr( ,N,LPr»KAf ,NtnAPF) 
*AMAS5 is RFPLACFP PY HCftF EHAPFS. 

♦♦STIF IS PtSTRPVFP. 

USF5 tlAFA,PCPMl»FlPNl ,1NVA. 

MAX N s SCO 

CALL MPPFIA (♦AMASS, ♦♦STIF ,W?.W, FFFCfF'tFPD. HAS, FUTAPF) 
♦AMASS IS PFPLACEIi FV MCr’F SHAPEF. 

♦*STIF IS rrSTRCVED. 

USE S BTAFAfDCPMl ,Eir.M tI^'VA. 

MAX M = SCO 

CALL MPDEIB ( ♦AMAS S , ♦♦FLFX , W? , W ,FRFC ,N , FPC ,KAF ,KUTAFE) 
♦AMASS IS PFPLACFO PV MPOE SHAPES. 

♦♦FLEX IS PbSTPrYFP. 

USES PTAPA,PCLM1,FIC;H1 ,IHV4,MLfLTA. 

MAX N = SCO 
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.13 PIGin Bcrv MCPfS 
.13.01 CALCULATION 

CALL RPTCl (XY2,XYZPFF,J»'Pf tJVFC.Z ,NNPnF,*NRZ,*NC2,KXJ, 
K2» 

♦NR2t^C^ WILL PE OEFINEr PY RE TGI . USE SYNEC L. 

USES RFVAOD. 

CALL RPTG2 I XRT ,X YZ RFF , JDOF, JVFC *7 tW^OOF ,^NRZ t *NCZ fKX J , 
K2I 

*NPZ,NCZ WILL PF OFF INFO PY RBTC?. USF SYMBOL. 

USES REV ADD. 

.13.02 PPTHC-NORMALIZAIION 

CALL CNPPM <PPMrOEfAMASS,N,NFPMrD,KPAI 
USES BTAB,F1CN1 ,MULTA. 

MAX N = 750 

MAX MRP MOD = 6 
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.14 INFPTIAL TRANSFORM AT I PM 

CALL UMAMl C AMASS »P ^MPOe ,Z ,N,NRPMi r ,KAFZ ) 

USFS PAPT*PTAP,INV1,MULTF 
MAX N = 250 

MAX M«BMCD = 6 
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.15 

INTFPNAL LOAOS 


.15.01 

PFAM - 5HFAR, 
CALL VMl 

HDMENT loads 

(XVEC ,DI57,CCNC,AMP1 , AMP?,tONVRT ,ZVECV,2VECH , 
NX,ND,NC,NA1 ,MA?,KO,kC,KA1 ,KA?» 

.15.0? 

PFAK - SHEAR, 
CALL VHTRl 
♦NZ=?NPP 

HOMENT TRANSFORMATION 
(PP,Z,NPP,*NZ,K2) 

HILL PE DFFINrr FY VMTRl . USE SYMEOL. 
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.16 time PFEPCNSE 

.16.01 FO*;rF IE riFTAIN'FD FY LlNFAP IMFFPCl AT 1ft.' USING TAFT.TAPF 

fALL TRSPl |aa,*F v*C ,*D, ■»#??, TAFF,XDOtXO»SIAPTT,PFLTAT» 

ENDT,NV^RITE ,NX ,NRTAFtNf TAB,6HXNAME tKAFCCt 
KTAF,NXTAPF,NrTl ) 

♦MATRICES A,F,f»P APE PFSIPCYFC'. 

USES INV1,MULTB. 

MAX NX = 250 
MAX NRTAF = 500 

CALL TRSPl E ( F ,f ,0 , T AF T, TAF F t XOC ,XC , ST ART! ,CF L TaT »E Nf T , 

NWP ITT ,NX ,NF T A F ,NC T AP , 6HXN AKE ,Kf f D »KT AE ,NXT APE ) 
MAX NX = 250 
MAX NPTAP = 5CC 

CALL TPSP2 (*A,^P ,vC,*r, TAFT, TAPF,XnC,>0, STAPH ,DFLIA1 , 

FMr<T,PFTA,NWPTTF,F'X *MRTAP,NCTAB,t>HXNAME ,KA'=^CDt 
KTAF, NXTAPE.NUTl) 

♦mat PICES A,P,C,D APF DFSTPUVFP. 

USES INV 1 , MULT, MU I TF. 

MAX NX = 25C 
MAX FS?TAP = 250 

CALL TPSP3 (♦AVFC,tPVFC,»CVFC,tr,TAPT,TAFF,^XDO,*XO, 
STAFTT,rELTAT,ENDT,NWP ITl , NX ,NRTAE , 

NCTAP ,6HXNAMiF ,Kn,KTAP ,NXTAPE ) 

♦VFCT(T?S, AVFC,PVrC,CVEC,XDC,XO, MATRIX D API DFSTROYEO. 

MAX NX = 250 
MAX NRTAF - 500 

MAXIMUM UMCUE TIMES PAST STARTT IN TABT = 250. 

.16.02 FCPCF IS CALCULATFO AS <1 - CPS 1/2 

CALL TPSPIA f^A,^F ,^C,^0, FMAG,FPP,VFL, CL, XOCt i'O, STARTT, 
PFLTAT,FNnT,NWRITf ,NX,NF,6HXNAMF ,KAFCD, 

NXTAPE ,NUT1 ) 

♦matrices A,P,C,D are DE^^TPrYFO. 

USFS INV1,MULTF. 

MAX NX = ?5C 
MAX NF = SCO 

CALL T»SP1C (F ,C ,L ,F MAC ,F F F ,Vr L , GL , XPO, XO, ST AR TT ,0F LT AT , 
ENfTTNWPITE ,NX ,NF ,6HXNAME ,KF CO, NXTAPE ) 

MAX NX = 250 
MAX NF = 5CG 

CALL TPSP2A ( *A ,*F ,*C ,tP , P MAG ,F P F ,VFl ,G L ,X00 ,XC, START T , 

OFL^AT ,F NOT, RE TA.,NWP IT F ,NX,NF,6HXNA ME ,KAfCf. , 
NXTAPE, NUTl) 

♦MATRICES A,f^c,n ARE nFSTRrvrr. 

USES INVl ,MUIT,MUITP. 

MAX NX = 250 
MAX NF = 250 
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.16 TIMF RFSPOMSF (CPNTINUFDI 

.16.03 ADDITIONAL F&UATICNf. 

CALL TRAF2 (6HXRUNNC,6HXK'AMF , I FA , A ,IF F ,P , JFC ,C , IFp, D, 

IFF ♦F,ZTMM,ST/PTT,EN0T*MLTXTP,NWPITt,*2IDFNT, 
♦♦STA,6HZNAMF ♦NZ,KZ,NXTAPF,N2TAPF , STORFZ > 

♦Z HEADING U2A6 FORMAT). 

♦♦STATIOF'S (A6 FORMAT). 

MAX NX = 250 

.16.04 TIMF RFEPCNSE MAX-MINS 

CALL TRMM (6HXRUMN0,6HXNAMF , XTMM , FT APTT ,FNLI , =»NX , 
KXfNXTAPF) 

♦NX IS DEFINED IN SUPRCOTINF, USE SYMPCL. 

MAX NX = 250 

.16.05 POWER SPECTRAL DENSITY OF ADEITIONAL EQUATIONS 

CALL TRPSD l6HXRUNNr,6HXMAMF , IR AE , lEXP t ST AFTT ,MLTXTP , 
ZPSntNFRFO,TlMPFR,NXTAPF,WRKV) 
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.17 fpfoueascy respowse 

.17.01 RESPONSE SOLUTICN 

CALL FRl f*Af*Bf*C,*0*TAPW,TAPF,lMEGAtNX,NOMEGA, 
NPTAP ,NCTAP tKAECn.KTAB *WPKM»NTAPE 1 
♦MATRICES A,e,C,n ARE OFSTPrYEr, 

P MUST PE NUN-SlNfl'LAR . 

USES INV1,MULT,MULTE. 

MAX NX = 50 

MAX NRTAE = 80 

.17.02 AOriTIONAL ECUATICNS 

CALL FPAEl (At*STA,^^21DFNT,NZ»NX»NrMFC AtKA.NXTAPF, 
NZIAPE) 

♦STATIONS IA6 FORMAT). 

♦♦Z HEADING ll2Afc FORMAT). 

MAX NZ = 80 

MAX NX = 50 
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APPENDIX P. rUMMAPY OF CALLING INSTPUCTirNS - 
FPARSF FORMA SOBPPOTINFS. 


♦♦♦ LIST OF SYMBOLS ♦♦♦ 

A = INMfT MATRIX 

ALFhA = INPUT SCALAR 
BETA = INPUT SCALAR 

NUTA = LOGICAL MIMBER OF UTILITY TAPE CONTAINING INPUT MATRIX A 

NUTB = LOGICAL NUMBER OF OTIIITY TAPE CONTAINING INPUT MATRIX B 

NUT2 = LOGICAL NUMPFP OF UTILITY TAPE CONTAINING OUTPUT MATRIX L 

AVFC = INPUT VECTOR (ROW OP COLUMN MATRIX) 

IVFC = INPUT VECTOR (ROW tR COLUMN MATRIX), INTEGER 

Z = RF5ULT MATRIX 

ZVFC = RESULT VECTOR 

W = MATRIX WORK SPACE 

NR = NUMBER OF ROWS 

NC = NUMBER OF COLUMNS 

KR = ROW 0 IMF n: ION 

KC = COLUMN DIMENSION 

V = VECTOR WORK SPACE 

LV = VECTOR WC-RK SPACE, INTFOFR 

KV = V,LV DIMFNSION 

NUTI = LOGICAL NUMBER OF ITH UTILITY TARf 

NTAPE - LOGICAL NUMBER OF FORMA RFS^RVE TAPE 



01 


MI«:CFLLANFril»S 


. 01.01 


, 01 . 0 ? 


,01.03 


,01. OA 


CONVFRSION 

CALL YSTOD (NUTA t Z t*NR .♦NC C ,V t LV , KV.NUTl ) 
♦NP.NC WILL BF DFFINEO BY YSTOO. USE SYMBOLS. 
USES Y1N,YIN1 ,Z2EPMb. 

OROFR MATRIX NCN-ZFRC ELEMENT LCFATIPNS ROWWISE 
CALL YLCiED (NUTAtVtLV,KV,NOTl ,NOT? ) 

USES Y1N,YINI,YCUT,YOUT1,YPART,Z2EOMB. 


ELIMINATE 
CALL ‘ 
US 


REPARTITION 
CALL VPART 


TE ZERO LEM ENTS 
VNrZfck (NUTA,V,LV,KV,NUT1 ) 

St 5 YIN ♦YINI,VnUT,YOUTI,YPAPT,ZZPnME, 

TT nm 


(NUTA»V»LVtKVfMUTI ) 


USE S YIN ,YINI ,YOUT,YOUTI ,Z2B0MB 


,01. OS COPY TO ANOTHER TAPE 

CALL FPUAL (NUTA,NUTZ tV»LV,KVI 
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•02 INPUT 

. 02,01 RFAL NlfMRFRS 

CALL YRFAP (NUTZ.V, LV,KVtNUTl I 

UFFS INT^PF,LTAPF,PAGFHn,YIN,VlNl,V(iUT,YClUTl,VPAPT, 
YRTAPF*YWRITF,YWTAPF,7ZP0MP . 

CARO INPUT I SANE AS RFAO FXCF'*T SHAPE) 

FIRST CAPO - 2NAMFfNP,NC WITH A6,14,l; FORMAT, 
MATRIX SHAPF IN COLUMNS 16-.' I . 

REMARKS IN COLUMNS 22-69, 

$ IN COL 7? FOR WRITF TAPF INITIALIZE, 
BLANK, RFWINO, LIST, OR TAPLID (FOR WRITE 

TAPE) IN COLUMNS 73-70, 

NWTAPF IN COLUMNS 79-60, 

MIDOLF CARDS - DATA WITh 2I5,AF17 FORMAT, 

LAST CARD - TFN ZEROES IN COLUMNS. I-IO, 

TAPF INPUT (SAME AS READ) 

ONE CARD - ZNAMF,C OP -LOCATION,NRTAPE (WITH - FOR 

NO WRITF OUT),ZRUNNO WITH 
A6,I5,A6 'OPMAT 

BLANK, REWIND, OR LIST (FOR READ TAPE) 

IN COLUMNS r2-27. 

REMARKS IN COLUMNS 2B-69, 

$ IN COL 72 FOR WRITE TAPE INITIALIZE, 
ELANK,PFWIND,L 1ST ,(?R TAPEID (FOP WRITE 

TAPF) IN COLUMNS 73-7P, 

NWTAPE IN C.fLl»MNS 79-PO, 
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.0? 

CUT Pin 


• 

O 

U) 

• 

o 

PRIMT RFAl. NUMeEPS 

CALL YWPITF (NUTA,6H ANAMF ,V,LV,KV) 
USES PAGFHDfYIMt /I^!I»Z?Pl’P'F . 


.03.02 

PUNCH PFAL NUMFEPS 

CALL YPUNCh INUTA,6H ANA>*F v </t LV , KV) 
US»T vin,y1N1,2ZBCMP. 
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.OA TAPE CPFRATIENS 

REWIND NPTAPF,NWTAPE AT PEGINNING DF MAIN PPDGPAM. 

SEE APPENDIX A FEE TAPE INITIALIZING, LISTING, AND UPDATING. 

.04.01 READING 

CALL VRTAPE l6HZRUN*Nr,6H 2NAPF,NUT2 ,V,LV,KV,*NRTAPE ,NUT1 ) 
♦NPTAPE MUST PE A REAO,WPITF TAPE ONLY. 

USES LTAPE ,Y1N,YIN1 ,YCUT,YPUT1,Z2FCMP. 

CALL YIN (NUTA,2VEf.,LS,LEI 

CALL YINI (NUT,1A,NS,NF| 

.04.02 WRITIW 

CALL VWTAPE IN'*TA,6H ANAME ,V,LV , KV,*NWT APE I 
♦WTAPE MUST PE A READ, WRITE TAPE ONLY. 

USES YIN ,VIN1 ,22PrMP. 

CALL YCIT <NfTA,AVEC,LS ,LFI 

CALL YDUTI <NUT, I A, NS ,NE ) 
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.C5 GFNERATICN 

.Cb.fll BASIC 

CALL YUNITY fNUTZf NP ,V,LV*KVI 
USFS YCUT.YPUTI. 

CALL YZFRO (NUTZ«NR,NCI 
USFS YCUT,YDUTI. 

.C5.C2 RAYLEIGH VECTORS 

CALL YPVl CNUT2,NtNU,V,LV,KV,Nirri ,^uT2fK'UT3l 

USFS YIN,YlNl,YLr»0,VC!UT,YOUTI,YPART,YTRANS,ZZBOMP, 
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.06 AlOFPPA 

.f6.ri fCAlAP pppouri 

CALL YAA lALPHA,NUTA,Nin?,V,LV,KV,NUTl,NUT?| 

U5FS YTNtYlNIfVLCRP,YNt ?FR,Yrirr,VrUTl,YPART,2?EC Mf . 

.06.0? ADDITirN, SUPTRArilCN 

CALL YAAPP tALPhA,NUTA,PFTA,NUTP,NU7ZfV,LVfKVfNUTl, 
NUT?) 

U5FS XL0RD»YIN,Y1N1 .YLCFO, YNP2FR ,YOUT,Yf'UTl tYPAFT, 
YSYMLH,Y<YMUH.Z2PrMF. 

.06.0? MULTlPLIf-ATION - BASIC 

CALL YMULT (NUT A ,NUTF ,NUT2 ,V .LV »KV ,NUT 1) 

USES YIN.YINI ,YLrRP,YNr2FR,Yri'T,Yt'llTl,YPART,YSYMLHt 
YSYMUH,Z2PrMP. 

. 06 . 0 A PL'LTIPLIC ATION - SPECIAL 

CALL YMULTl (NUTA ,NUTE ,NU72 tV ,LV ,KV ,NUTI » 

USES YIN,YIN1 ,YL0Rr,YNU2FR,YrUT,YrUTI,YPART,YSYMLH, 
ySYMUHtZZFOME . 

CALL YMULT? ( M»TA »MUTF ♦NUT2 ,WI ,W2 , V ♦ I V , K V,KW ,NUT 1 ) 

USES YnTC!S»YIN,YlNI,YLrPr,YNn2FP ,YPUT,YOUT 1 .YPAF T ♦ 
YSYPLUtYSYMUHtZZFrME. 

CALL YML1T4 (NliTA ,E ,NUT2 .W,V,LV,KV,KWfNUTn 

USES YIN.YIM *YLCRC,Yf UT,YrUTI*YPART»Z2PCME . 

.06.05 TRIFLE MATRIX RFfDl’CT 

CALL YBTAE (NU7A,NUTE,NUT2,V,IV,KV,NUT1,NUT2> 

USES YIN tYINI ,Y LORD, YMULT tYNPZER ,YTUT ,YCU f I ,Y FART, 
YSYMLH,YSYMUH,2ZRrMF . 

.06.06 nFCCMPCSITirN 

CALL YDCM3A INUTA,NUTU,NUTr ,V,LV ,KV , NUTT , NUT? ) 

USES YIN ,YINI ,YLCPr,YmiT,YruTI,YPART,Y7P AMS ,2ZPf ME. 
CALL YDC0M3 (NUTA,NU7U,NUTn ,V ,l V ,K V ,NUT 1 ,NUT? I 

USES YIN ,YIM,YLrRD,YDUT,YOUTI,YPART,YTRANS tZZPCMP. 

.06.07 PACK S0LU7ICN 

CALL YFSL3A (MUTU,MU7r',NU7p ,NUT2 ,V, LV, KV,NUT1 ,NUT? ) 

USE S YlN,YINI,YLrRD,YrUT,YrUTI,YpART,YTRANS. 
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.07 POPIFICaTIPN 
.07.01 ASSE^'FLY 

CALL YASSPH INUTAtlPZ.JCZ , VtLV,KV,MllTl I 

♦FF 5UPF Z IS PFFINFO. (CALL YZFPC TP CLFAF NUTZ). 

•JSFS YIN,YI»iI,YLCRC,YCUT,YC’UTl,YPAKT ,YSYPLF t VSYPUH* 
ZZPPMF . 

.C7.C? DISASSFHFLY 

CALL YDISA (NL'T A « IP A. JC A ,Nl'TZ ,NF Z t NCZ , V,LV,F V,Nl'T 1 I 
USES YIN,YIM ,YLC:RD,VrUT*YCT!Tl ,VFART,Z2EPHE. 

.07.0? REVI5ICF’ 

CALI YRFVAD ( ALPHA »r<UTA, I VFC , JVFC ,*NLI7Z ,V ,LV ,KV ,NlTl , 
NUT?,N't'T3,NUTi:.) 

*PF Sl'RF 2 IS OFF INFO. (CALL YZFPP TC CLEAR Nl'TZI. 

I'f FS XLPFn,YAAFP,YIN,YIMI,YLnPr',YNn2FP,YnUT,YnUTl , 

YPAPT ,YSYPLH,YSYMl*H,2ZFrE^E. 

CALL YRVADl ( Al PHA , A , IJVEC ♦♦NPTZ »NP A,V, LV t KV,KA .NMT 1 ♦ 
NI'T?,NIjT?,NUT4) 

♦PE SIPF 2 IS OFFINED. (CALL YZFPP TC CLEAR NUTZI. 

Cf F S XLCP r »YA AE F , Y1 N , Y INI , Yl CFC, YNP2 f F , YPL T , YPCT I , 
YPART,Y5YPLH,Y5YMIiH,Z?FrPiP, 

CALL YRVAD? ( NIT A ,MiT2 ,NP Z ♦ W,KW , V»L V tK V ,NUT 1 1 NUT2 tH'l ? ) 

U5F5 YIN.YINI , YCUT, YCHT 1 ,YPAP T,Z ZPC'MF . 

CALL Vf VAD3 (NCTA.NCTZ ,NP? *NC2,W ,KW,V,LV,KV,NUT1 ,Nl<T?,NUT?l 
USES YIN,Y1N1 ,YCT!T,Ya'Tl,YFAPT,/ZFrP!E. 

CALL YRVISI (A.JVrC.NUTZ.NFAZtNCAtNCZ^V.LV.FV.KAj 
U5 F S XL 0 P D f VOLT , Y f UTI » 2 7 FCMP . 

CALL YPVTCD (NIT A. IVEC ♦ JVEC ,Z ,NR Z ,NC2 ,V ,LV,KV,K2 I 
USES YIN ,Y1NI fZZFC'HB. 

.07.0A TPAMSPCSE 

CALL YTFUVNS (NUTA ,MIT2 t V ,LV ,KV tM'TI tM'T? » 

USES YIN , YINI tYLCRD ,YCl*T ,YrUT I ,YPART ,Z2FrPF . 

.07.0* SYPMFTPIZF 

CALL YSYNLEi (NUTA2 ,V,LV,KV,NIJT1 ,NLiT?) 

USES YIN, YINI ,YLnFn,YNr’2EP,YPlfT,YnL’T I,YPAP1 . 

CALL YSYMUH (NUTaZ ,V,LV,KV,NI»T1 ,N(fT? I 

USES YIN, YINI ,YLrRr>,YNr'ZFR,yruT,Yrt!Tl,YPAPT. 

.C7.0A NULLIFY 

CALL YZFPLH (NUTAZ ,V,LV,KV,NUT1 , NUT?) 

U5FS YIN, YINI ,YLCPP,YrttT,YrUTT,YPART ,221 ( ME . 

CALL YZFFUH ( MH AZ ,V ,LV,K V,NUTI , NUT ? ) 

USE 5 VIEf ,Y1NI ,YirPr,Yf UT ,YfUT I ,YF'AeT .ZZEE NF . 

,07,07 DIAnCNALlZE 

CALL YDIAO ( NUT A,NUT2 , V , LV ,KV) 

USES YIN, YINI, YCUT, YCUTT,7ZPf'PF, 
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• CP 5.7TFFFJFSS 
.rp.Cl PEmiCTION 

CALL YSPFD2 (NUT A,NU .K ,NUT7 ,NP ♦ I FI , V ,LV »KV »NU1 1 ♦NL'7 2 ♦ 
NUT3,NUT4) 

USFJ: Ya<;SFN,VD 1SA,Y1N,V1NI ♦YLnF'>,YNC2FP,YOUT,Yrini, 
YPAFT,YSYNLH,Y!;YNl!H,YTPAN5tYUNITY,Y2FRr,2ZPrMP 
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.r^ MODAL 

.09.01 itfpative payleich-fitz 

CALL VMfOPA (N'UTM,NUTK,NUTZtW?,W,FPrO,NW,NUtSHIFT,TrLZ t 
TDLW?,MAXIT,lFPRMT,V,LV,AtS,KVIN,KA, 

NUTl ,NUT2 .NUT3 ,NUT4 ,NUT!>,NlfT6 I 
USFS PTABA2,EIGNlA,INV4,M0rElX,NAME.PAGEHD,TIMCEK, 
WFITF,WRITTM,XLORD,YAAEE,YE«L3A,YOCM3A,YOTCJ ,YIMf 
Y IM I ,YLPR D , YMUL T1 , YMUL T? , YMUL T4 , YN( Z EP , YOUT , Yt UT I , 
YPART,YR Vl,YS10DtYSYMLH,YSYMUH,YTRANS,YWP nE,ZZPPMP. 



APPFNPIX C. rP CALL INC INSTPICT IONS - 

FINITE ELEMENT RCDTINFS. 


A CH ARfUMFNT MAY ALSC PE A VARIAPLE READ WITH AN A6 FORMAT CF 
rETAINED WITH A DATA STATEMENT. 


♦♦♦ LIST GF SYM6CLS »♦» 


NUTE 

NUTFL 

NUTK 

NUTKX 

NUT IT 
NUTM 
NUTMX 
Nl»TST 


LGGICAL NL!MFER CF UTILITY TAPE ON WHICH ELEMENT UNIT LOAO 
PUCKLINf MATRICES AND TVECS ARF GirTPia. 

Nl»TE MAY PE ZEPG IE EUCKLINC MATFICFS ARE Nt»T FORMED. 

USES FORTRAN READt WRITE. 

LPCICAL NUMPER OF TAPE CONTAINING EHMENT INPUT DATA FOP 
THIS SUFFOUTINE AND SUPROUTINES AXiALt ETC GIVEN FY NAMEL. 
IF NUTEL = 5f DATA WILL PE READ FROM CARDS. 

LOGICAL NLIMPFR OF UTILITY TAPE ON WHICH ASSEMELEO 
STIFFNESS MATRIX IS OUTPUT IN SPARSE NOTATION. 

NUTK MAY FE ZERO IF STIFFNESS MATRIX IS WT FORMED. 

USES FOPMA YINt YOUT. 

LOGICAL NUMBER OF UTILITY TAPE ON WH ICH ELEMENT 

stiffness matrices (Same as clopal loads transformation 

MATflCtS) AND IVFCS A«^r STCRFD. 

NUTKX MAY FE ZERO IE STIFFNESS MATRIX IS NOT FORMED. 

USES FOP TRAN READ. WRITE. 

logical NUNPER PE UTILITY TAPE ON WHICH TLEMENT LOCAL 
LOAD transformation MATRICf^ AFfO IVECS APE OUTPUT. 

NUTLT MAY'VE ZERO IE LOAD T*- , ' SEOPMA TTONS ARE NOT FORMED. 
USES FORTRAN READ, WRITE. 

LOGICAL NUMBER OF UTILITY TAPE ON WHICH ASSEMELEO 
MASS matrix is OUTPUT IN SPARSE NOTATICN. 

NUTM MAY PF ZERO IF MASS MATRIX IS NOT FORMED^ 

USES FORMA YIN, YfUT. 

LEGICAL MUMPER OF UTILITY TAPE ON WHJCH FLFMFNT 
MASS MATRICES ANO IVECS ARE STORED. 

NUTMX MAY PE ZERO IF MASS MATRIX IS NOT FORMED. 

I’SES FORTRAN RFAP, WRITF. 

LOGICAL M'MPER OF UTILITY TAPE ON WHICH ELEMENT 
STRESS transformation MATRICOS IVECS APE CUTFi:T. 

NUTST MAY BE ZERO IE STRESS TR ANS FfR MAT I rtvS ARE NO'T FORMED. 
USES FORTRAN REAP, WRITE. 



Cl 


MISCFLLANFOUS 


.Ol.Cl DIFFCTIDN COSINES 

CALL DCOSIA (CJtFJ.ZtKCJtKFJfKZ ) 
USFS FULEP,MULTF,ZZPCMP. 

CALL DCnSlP aJtFJ,Z,KCJ,KFJ,KZI 
USES FULEF.MULTP.ZZPOMP, 

CALL DC052 (C J* FJ,Z ,KC J,KF J,KZ ) 
USFS ElJI.FR,MULTP,ZZPrMP. 

CALL DC0S3C (CJtFJfZtKCJ«KEJtKZ ) 
USES EULEPtMULTPfZZPOMP. 

*01. C? FUIFP ANCLES 

CALL EULFR fF,P,KR| 

• 01.03 TETPAHEPRCN GFOMETRY 

CALL TEGEON 1C J, JM ,VL ♦OV,KC J, IFB AD) 
USES vcROSStVnrT. 
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.02 FINITF ELEMENT MASS, S'.IFFNFSS, PUCKLINO, LOAD TP ANSFOPMATIPN, 

STPFSS TPANSFCPMATIPN 

CALL FINFL (> YZ , JOCF ,FUL ,NUTFL ,NJ,NUTM ,NUTK ,NUTl T, 

NUTST,NIITP,V,LV,KV,KRX ,KR J ,KP.F ,NUTMX,NUTKX, 
NUTl ,NUT2,NUT3I 

USES AXI At, FAR, FLUID, GPAVTY, PACE HP, QUAD, PFCTSP,TPNCL, 
YPVAD2,22Fr!ME. 

CALL AXIAL ( XYZ , JDPF ,EUL ,NUTEL , NJ, NUTMX ,NUTKX ,NUTLT, 
NUTST,W,T,S,KX,KJ,KE,KW| 

USES ATXPAl ,ATXPP ,0C0S1 A ,EUL E P ,K T Al ,M1A1 ,MlA2,MASl A, 
MULTA,MUlTP,PAGEMn,STFlA,Z2PnMB. 

CALL BAP IXYZ,J0PF,FUL,NUTFL,NJ,NUTMX,NUTKX,NUTBX, 

NUTl T ,NUTST , W , I,S,KX,KJ,KE,KW» 

USFS ATXPA1,P1AI,B1A2,PTABA,BUCIP,0CDS1P,EULER,KIa1, 

KIP 1 ,K1C1,MJ A1 , Ml A?, Ml PI ,M1P?,M1C1 ,M1C?,MAS1B, 
MULTA,MULTP,PAGFHP,STF1E ,ZZEnME . 

CALL FLUID IX YZ • JOOF , FUL ,NUTFL , NJ,NUTMX,N«m X,NU7LT,NUTST , 
W,T,S,KX,KJ,KE,KW) 

USES TEGEOM,VCROSS,VDOT,ZZPrMP. 

CALL GF AVTY ( X YZ , JD^JF, EUL ,NlrTC L , NJ,MUTK X, W, 1 , S ,KX,K J,KE,K W ) 
USFS KGRAV,MULTA,MULTB,VCRrSS ,ZZEOMB. 

CALL OUAD fXYZ,JOnF,EUL ,NUTfL,MJ,NUTMX,NUTKX,NUTBX, 

NUTLT ,NU1ST,W,7,S,KX,K J,KE,KW) 

USES ATXEAl,PTABA,DCOS2,niLFR,K2Al,K?El,M2Al, 

M?A2,M2P1 ,M2B2,MAS?,MAS3,MULTA,MUL7F,PAGFHD, 

REV ADP,S7F2,S7F3,ZZ POMP, 

CALL PECTSP (XYZ,JPnF ,EUL,MUTEL , NJ,NUTMX ,f'UTKX , NUTLT, NUTS7 , 
W,T,S ,KX,KJ,KF,KW) 

USES M A 5 3 A , PAGEHD , ’^TF ? A , ZZ POM F . 

( *.LL 7RNGL ( XYZ , JDOF ,EUL ,NUTEl , NT,Ni:TMX ,NUTKX ,NU7P X, 

NUTLT, NIITST,W,T,S,KX,KJ,KF ,KWI 
USFS ATXP AI ,ETAPAtDC0S2,FULEP ,K?Al ,K2BI ,M2Al , 

M2A2,M2P J ,M2f 2,MAS2,MULTA,MULTP ,PAGFH0,STF2 ,ZZPUMP . 
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.03 MASS 
.03.01 por 

CALL MASlA I C J. FJ ,A1 , A2 * Rr,6MMA KFM tKCJ»KEJ*K2) 

USFS AlXPP,EULEPfMlAl,MlA2,Z2POMP. 

CALL MlAl (Al,A2»RL,Pr,2,KZ) 

CALL MIA? (Al,A?,PL,PO,2,K2) 

CALL MlCl fPIl.PI2tRL»Rr,Z,KZ) 

CALL MIC? IPIl tP12,RL»PP,7,K2) 

.03.0? BEAM 

CALL MASIP (CJtFJ,Al,A2»PIl,PI?,Kr,6HNAMEM t 2 ,W »KC J,KF J, 
KZtKWI 

USES PTAPA,DCPS1B,FULER,M1 Al , m7a?,M 1E1,M1B2,M1C1 , 
M1C?*MULTP ,Z2PnHP . 

CALL MlPl Ul,A?,RL,Rr,2,K2 ) 

CALL M1F2 <A1,A?,RL,RC,2,K2) 

.03.C3 TRIANGLE 

CALL MAS2 (CJ»FJ,TMAS ,Rf .6HNAMEM fZ,Wl ,W2,KC J,KE J,K2 ,KWl , 

KW? ) 

USFS FTAE AtDCrS?. FULFP ,M?A1 ,M2A? .M2P 1 tMZB? »MULTP »2 ZFCMP . 
CALL M2AI IX2.Y3 f7H,R0tZ*K2 ) 

CALL M2A? (X2,X3»Y3,THtRnt2tT,R,KZ.KT,KRI 
USES PTAPA. 

CALL M2F1 lX2tY3,TH,Rr,2*KZ) 

CALL M2P2 (X2»X3,Y3,TH,Pr,Z,T,R,KZ,KT,KP.) 

USFS PTAPA. 

.03.04 OUAPRILATEFAL 

CALL MASS <CJ,FJ,TMAS »R( tCHNAMPM .Z»W1 ,W2,KCJ,KEJ,K7 .KKl , 

KW2 ) 

USES PTAPA,DCrS2TFUlER,M?Al,M2A2,M?Pl,M2P2fMAS2, 
MULTP,Z2f OMP . 

,03. OS PFCTANfULAR SHEAR PANEL 

CALL MAS3A (C J , E J t TMAS ,Rr, fcHMAMEM , 2, W1 , W? ,KC J,KEJ , 

KZ.KWl ,KW?) 

USES E.TAPA,DCrS 3 C»M 3 Cl ,Z?PPMP , 

CALL K3fl (X3tY3,TH,P0tZ,KZ J 





.04 STIFFNFS5 
.04.01 ROD 

CALL STFIA ( C Jt F J ?A 1 , A? , f t^HNAMFK ,6HF' AMF S T , S , 7 L fTS.NF ST , 
KCJfKFJ ♦KS,KTL,KTS » 

USFS ATXE Al,DCrSlA,FULFF ,K1A1 ,MUL7A ♦ HULl F .ZZFOMP . 

CALL KlAl (AltA?.PL»FtZfT5,K2,KT5 ) 

CALL Kiel (TJl,TJ2tPltP?*'!L,G,7,TS*KZ,KTS1 

.pA.r? PFAM 

CALL STFIP (rj,FJ,Knrt,Al*A2»TJl,TJ?,PI21*eiZ2*riYl, 

PIY?,R1 ,P?,CV1,CY2,C21 ,C22,SF,F,f ,6HNAMFK , 
bHNAMFST,S,TL,Tf.fNPST,KCJ»KFJ,KS,KTL,KTS) 

USFS aTXBAI ,DCrSlP,KlAl,KlPl ♦KIC 1 ♦MUITA,F»UL7P t ZZHOMP. 

CALL KlPl (PI1,PI2,C1,CI*A1,A2»SF,RL,F,G,Z,TS,K2,KTSI 

.04.03 TRIANGLE 

CALL S7F2 tC Jt F J ,TRlFM , TFFNt F . AMJ, GM^AMEK ,6hNAMEST,S ,TL ,7 S ♦ 
NPST,KCJ,KFJ,KS,KTL ,K7 S ) 

USES ATXPAl ♦rcrS2,FULEP,K2Al ,K2P1 t^^^»L7A»MULTP^2ZPnMP. 

CALL K2A1 (X2tX3,Y3,TH,f ,AMU.Z,T,R,KZ,KT*KR) 

USFS PTAPAtMUL7A,Z2POMP. 

CALL K2F1 (X2,X3,Y3»TH,F,ANU,2 ,TS ,T,KZ,KTS,KT) 

USFS PTAPAtMULTA,Z2FCMF. 

.04.04 OUADPILATFrAL 

CALL STF3 (CJ,E JtTMEM.TPfN.F, ANU,4HNAMFK ,6HNAMFST,S ,TL t7S ♦ 
NFST,KCJ,KEJ,KS,K7L ,K7S) 

USES ATXPAl ,DCUS2tEULFF ,K?A1 , K2P 1 »MULTA , MULTE ,RFVAr>D, 
STF2,ZZPDMP. 

.04. OS PFCTANGULAfi SHEAR PANEL 

CALL STF3A ( C J t F J tTH, f *6KMMEK .6HNAMF ST, S ♦TLfTSfNRST 
KEJ,KS,KTL,K7S I 

USFS ATXPAl ,DCPS3C,K3C1 ♦MULTa.ZZPUMP. 

CALL K3C1 (X3,Y3,TH,G,Z,7,K7,KT) 

.04.06 FLUID 

CALL KGPAV (C J , JM ,F V , A ,W ,KW,KC J I 
USES MULTA,MULTe,VCRCSS. 
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.C5 PUCK LING 


CALL PUCIP (CJtFJ,KrDEB,6HNAHEB *Z ,W,KCJ,KE J,KZ *KW » 
USE S FI A1 »FlA2fP.TAPAtOCnSlP,EULEP tMULTB.ZZhdMP. 
CALL FlAl <PLtZtKZ) 

CALL B1A2 (PLtZ.KZ ) 
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.06 FLUID PRFSSURF 

CALL PRFSS (C Jf T,NJD,NCrL,KCJ,KW) 
USF5 PFVADDtlNVlNP.MltLTP. 



